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Abstract. Gamma-Ray Burst (GRB) explosions from the first generation of stars 
offer an exciting opportunity to probe the epoch of reionization. Clues about how and 
when the intergalactic medium was ionized can be read off their UV emission spectrum. 
Roughly a percent of all GRBs should be strongly gravitationally lensed by intervening 
stars. A microlensed lightcurve can be inverted to reconstruct the surface brightness 
profile of the GRB image on the sky, with micro-arcsecond resolution. 



1 Introduction 

Since their discovery four decades ago, quasars have been used as powerful light- 
houses which probe the intervening universe out to high redshifts, z ~ 6 (|25[p[). 
The spectra of almost all quasars show strong emission lines of metals, indicating 
super-solar enrichment of the emitting gas |24j . This implies that at least in the 
cores of galaxies, formation of massive stars and their evolution to supernovae 
preceded the observed quasar activity. If Gamma-Ray Bursts (GRBs) originate 
from the remnants of massive stars (such as neutron stars or black holes), as 
seems likely based on recent estimates of their energy output f45| , [iT {lO | , then 



they should exist at least out to the same redshift as quasars. Although GRBs are 
transient events, their peak optical- UV flux can be as bright as that of quasars. 
Hence, GRBs promise to be as useful as quasars in probing the high-redshift 
universe. 

Not much is known observationally about the universe in the redshift interval 
z = 6-30, when the first generation of galaxies condensed out of the primordial 
gas left over from the Big Bang (see reviews || and j30|). Observations of the 
microwave background anisotropics indicate that the cosmic gas became neutral 
at z ~ 1000 and remained so at least down to z ~ 30 (see, e.g. |j43f). On the 
other hand, the existence of transmitted flux shortward of the Lya resonance 
in the spectrum of the highest-redshift quasars and galaxies (see, for example, 
Figure , indicates that the intergalactic medium was reionized to a level better 
than 99.9999% by a redshift z ~ 6. This follows from the fact that the Lya 
optical depth of the intergalactic medium at high-redshifts (z ^> 1) is 



TV, 
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where H « lOO/i km s 1 Mpc 1 fiH 2 (l + z) 3 / 2 is the Hubble parameter at the 
source redshift z, f a = 0.4162 and \ a = 1216A are the oscillator strength and 
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the wavelength of the Lya transition; uhi{z) is the average intergalactic density 
of neutral hydrogen at the source redshift (assuming primordial abundances); 
fi m and J7f, are the present-day density parameters of all matter and of baryons, 
respectively; and xhi is the average fraction of neutral hydrogen. Modeling Q 
of the transmitted flux in Figure [l] implies r Q < 0.5 or xhi < 10~ 6 , i.e., the 
low-density gas throughout the universe is fully ionized at z ~ 6! 
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Fig. 1. Optical spectrum of a z — 5.8 quasar, discovered by the Sloan Digital Sky 
Survey (Fan et al. 2000). 



The question: how and when was the universe reionized? defines a new fron- 
tier in observational cosmology |M| . The UV spectrum of GRB afterglows can be 
used to probe the ionization and thermal state of the intergalactic gas during the 
epoch of reionization, at redshifts z ~ 7-10 [Q. The stretching of the temporal 
evolution of GRB lightcurves by the cosmological redshift factor (1 + z), makes 
it easier for an observer to react in time and take a spectrum of their optical-UV 
emission at its peak. 

Energy arguments suggest that reionization resulted from photoionization 
rather than from collisional ionization 1^,0 . The corresponding sources of UV 
photons were either stars or quasars. Recent simulations of the first generation 
of stars that formed out of the primordial metal-free gas indicate that these 
stars were likely massive If GRBs result from compact stellar remnants, 

such as black holes or neutron stars, then the fraction of all stars that lead 
to GRBs may have been higher at early cosmic times. This, however, is true 
only if the GRB phenomena is triggered on a time scale much shorter than 
the age of the universe at the corresponding redshift, which for z 3> 1 is ~ 
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5.4 x 10 s yr (h/Q.l)- 1 (Q m /Q.i)- 1 l 2 [{l + z)/10]~ 3 / 2 . This condition may not 
hold, for example, for neutron star binaries with an excessively long coalescence 
time. 

In §2 we briefly discuss the expected properties of GRB afterglows at high 
redshift. We then describe the use of distant GRBs for two studies: (i) probing 
the intergalactic medium (IGM) during the epoch of reionization (§3); and (ii) 
finding intervening stars at cosmological distances through their gravitational 
lensing effect (§4). 



2 Properties of High-Redshift GRB Afterglows 

Young (days to weeks old) GRBs outshine their host galaxies in the optical 
regime. In the standard hierarchical model of galaxy formation, the character- 
istic mass and hence optical luminosity of galaxies and quasars declines with 



increasing redshift [p3|j41| ,|3[| . Hence, GRBs should become easier to observe than 
galaxies or quasars at increasing redshift. Similarly to quasars, GRB afterglows 
possess broad-band spectra which extend into the rest-frame UV and can probe 
the ionization state and metallicity of the IGM out to the epoch when it was 



reionized at a redshift z ~ 7-10 1 30 . Lamb & Reichart []29j have extrapolated 
the observed 7-ray and afterglow spectra of known GRBs to high redshifts and 
emphasized the important role that their detection might play in probing the 
IGM. Simple scaling of the long-wavelength spectra and temporal evolution of 
afterglows with redshift implies that at a fixed time lag after the GRB trigger 
in the observer's frame, there is only a mild change in the observed flux at in- 
frared or radio wavelengths as the GRB redshift increases. Ciardi & Loeb [Q 
demonstrated this behaviour using a detailed extrapolation of the GRB fireball 
solution into the non-relativistic regime (see the 2/im curves in Figure ^|). De- 
spite the strong increase of the luminosity distance with redshift, the observed 
flux for a given observed age is almost independent of redshift in part because 
of the special spectrum of GRB afterglows (see Figure ||) , but mainly because 
afterglows are brighter at earlier times and a given observed time refers to an 
earlier intrinsic time in the source frame as the source redshift increases. The 
mild dependence of the long-wavelength (A b s > 1/im) flux on redshift differes 
from other high-redshift sources such as galaxies or quasars, which fade rapidly 
with increasing redshift p^ , |il] ,p[ . Hence, GRBs provide exceptional lighthouses 
for probing the universe at z = 6-30, at the epoch when the first stars had 
formed. 

Assuming that the GRB rate is proportional to the star formation rate and 
that the characteristic energy output of GRBs is ~ 10 52 ergs, Ciardi & Loeb 
predicted that there are always ~ 15 GRBs from redshifts z > 5 across the 
sky which are brighter than ~ 100 nJy at an observed wavelength of ~ 2/im. 
The infrared spectrum of these sources could be taken with future telescopes 
such as the Next Generation Space Telescope (planned for launch in 2009; see 



littp:/ /ngst.gsfc.nasa.gov/| ) , as a follow-up on their early X-ray localization with 



the Swift satellite (planned for launch in 2003; see http://swift.sonoma.edu/) 
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Fig. 2. Theoretical expectation for the observed afterglow flux of a GRB as a function 
of its redshift (from Ciardi & Loeb 2000) . The curves refer to an observed wavelength of 
5000 A(thin lines) and 2fim (thick lines). Different line types refer to different observed 
times after the GRB trigger, namely 1 hour (solid line), 1 day (dotted) and 10 days 
(dashed). The 5000A flux is strongly absorbed at z > 4.5 by intergalactic hydrogen. 
However, at infrared and radio wavelengths the observed afterglow flux shows only a 
mild dependence on the source redshift. 



The redshifts of GRB afterglows can be estimated photometrically from ei- 
ther the Lyman limit or Lya troughs in their spectra. At low redshifts, the 
question of whether the Lyman limit or Lya trough interpretation applies de- 
pends on the absorption properties of the host galaxy. If the GRB originates 
from within the disk of a star-forming galaxy, then the afterglow spectrum will 
likely show a damped Lya trough. At z > 6 the Lya trough would inevitably 
exist since the intergalactic Lya opacity is > 90% (see Figure 13 in [ftl[). In- 
terestingly, an absorption feature in the afterglow spectrum which is due to the 
neutral hydrogen within a molecular cloud or the disk of the host galaxy, is likely 
to be time-dependent due to the ionization caused by the UV illumination of 
the afterglow itself along the line-of-sight |TjJ • 

So far, there have been two claims for high-redshift GRBs. Fruchter [ |l2[ 
argued that the photometry of GRB 980329 is consistent with a Lya trough 
due to IGM absorption at z ~ 5. Anderson et al. Q inferred a redshift of 
z = 4.5 for GRB 000131 based on a crude optical spectrum that was taken 
by the VLT a few days after the GRB trigger. These cases emphasize the need 
for a coordinated observing program that will alert 10-meter class telescopes to 
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take a spectrum of an afterglow about a day after the GRB trigger, based on a 
photometric assessment (obtained with a smaller telecope using the Lyman limit 
or Lya troughs ) that the GRB may originate at a high redshift. 

In the following two sections, we use two examples to illustrate the usefulness 
of GRB afterglows for cosmological studies. 




Fig. 3. Sketch of the expected spectrum of a source at a redshift z s slightly above 
the reionization redshift z T eion- The transmitted flux due to HII bubbles in the pre- 
reionization era, and the Lya forest in the post-reionization era, are exaggerated for 
illustration. 



3 Probing the Intergalactic Medium 

The UV spectra of GRB afterglows can be used to measure the evolution of 
the neutral intergalactic hydrogen with redshift. Figure || illustrates schemat- 
ically the expected absorption just beyond the reionization redshift. Resonant 
scattering suppresses the spectrum at all wavelengths corresponding to the Lya 
resonance prior to reionization. Since the Lya cross-section is very large, any 
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transmitted flux prior to reionization reflects a large volume of ionized hydrogen 
along the line-of-sight. If the GRB is located at a redshift larger by > 18% than 
the reionization redshift, then the Lya and the Ly/3 troughs will overlap. Unlike 
quasars, GRBs do not ionize the IGM around them; their limited energy supply 
~ 10 52 ergs p5|pdpp| can ionize only ~4x 10 5 M Q of neutral hydrogen within 
their host galaxy. 

Quasar spectra indicate the existence of an IGM metallicity which is a frac- 
tion of a percent of the solar value || . The metals were likely dispersed into the 
IGM through outflows from galaxies, driven by either supernova or quasar winds 
(HJOl- Detection of metal absorption lines in the spectrum of GRB afterglows, 
produced either in the IGM or the host galaxy of the GRB, can help unravel 
the evolution of the IGM metallicity with redshift and its link to the evolution 
of galaxies. 



4 Cosmological Microlensing of Gamma-Ray Bursts 

Loeb & Perna |u| noted the coincidence between the angular size of a solar-mass 
lens at a cosmological distance and the micro-arcsecond size of the image of a 
GRB afterglow. They therefore suggested that microlensing by stars can be used 
to resolve the photospheres of GRB fireballs at cosmological distances. (Alterna- 
tive methods, such as radio scintillations, only provide a constraint on the radio 
afterglow image size [ p"8p5| ] but do not reveal its detailed surface brightness dis- 
tribution, because of uncertainties in the scattering properties of the Galactic 
interstellar medium.) 

The fireball of a GRB afterglow is predicted to appear on the sky as a ring (in 
the optical band) or a disk (at low radio frequencies) which expands laterally at a 
superluminal speed, ~ rc, where r 3> 1 is the Lorentz factor of the relativistic 



blast wave which emits the afterglow radiation [j44 39 3J3|,[20|]. For a spherical 
explosion into a constant density medium (such as the interstellar medium), the 
physical radius of the afterglow image is of order the fireball radius over r, or 
more precisely 12^] 



R s = 3.9 x 10 16 i^j H ° [jf^ 5/S cm, (2) 

where -E52 is the hydrodynamic energy output of the GRB explosion in units of 
10 52 ergs, ni is the ambient medium density in units of 1 enr 3 , and idays is the 
observed time in days. At a cosmological redshift z, this radius of the GRB image 
occupies an angle 8 S ~ R s /Da, where Da{z) is the angular diameter distance 
at the GRB redshift, z. For the typical cosmological distance, Da ~ 10 28 cm, 
the angular size is of order a micro-arcsecond (/ias). Coincidentally, this image 
size is comparable to the Einstein angle of a solar mass lens at a cosmological 
distance, 



/4GM^\ 1/2 ( M lens \ 1/2 ( D y 1/2 
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Fig. 4. A typical broken power-law spectrum of a GRB afterglow at a redshift z = 1 
(from Granot & Loeb 2001). The observed flux density, F v , as a function of frequency, 
v, is shown by the boldface solid line at an observed time tdays = 1 for an explosion 
with a total energy output of 10 52 ergs in a uniform interstellar medium (k = 0) with 
a hydrogen density of 1 cm -3 , and post-shock energy fractions in accelerated electrons 
and magnetic field of e c = 0.1 and es = 0.03, respectively. The thin solid line shows the 
same spectrum when it is microlensed by an intervening star with an impact parameter 
equal to the Einstein angle and Ro = [0 S (1 day)/#E] = 1. The insert shows the excess 
evolution of the break frequencies v\> = v a , v m and v c (normalized by their unlensed 
values) due to microlensing. 

where Mi cns is the lens mass, and D = D 0S D \/ D\ s is the ratio of the angular- 
diameter distances between the observer and the source, the observer and the 
lens, and the lens and the source pQ ]. Loeb & Perna (1998) showed that because 
the ring expands laterally faster than the speed of light, the duration of the 
microlensing event is only a few days rather than tens of years, as is the case 
for more typical astrophysical sources which move at a few hundred km s" 1 or 
- lCT 3 c. 

The microlensing lightcurve goes through three phases: (i) constant magni- 
fication at early times, when the source is much smaller than the source-lens 
angular separation; (ii) peak magnification when the ring-like image of the GRB 
first intersects the lens center on the sky; and (iii) fading magnification as the 
source expands to larger radii. 
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Fig. 5. The surface brightness, normalized by its average value, as a function of the 
normalized radius, r, from the center of a GRB afterglow image (where r = at the 
center and r = 1 at the outer edge). The image profile changes considerably between 
different power-law segments of the afterglow spectrum, F v oc v 13 (see Figure ^). There 
is also a strong dependence on the power-law index of the radial density profile of the 
external medium around the source, p oc R~ k (taken from Granot & Loeb 2001). 



Granot & Loeb jl9) calculated the radial surface brightness profile (SBP) of 
the image of a Gamma-Ray-Burst (GRB) afterglow as a function of frequency 
and ambient medium properties, and inferred the corresponding magnification 
lightcurves due to microlensing by an intervening star. The afterglow spectrum 
consists of several power-law segments separated by breaks, as illustrated by Fig- 
ure H The image profile changes considerably across each of the spectral breaks, 
as shown in Figure [5]. It also depends on the power-law index, k, of the radial 
density profile of the ambient medium into which the GRB fireball propagates. 
Gaudi & Loeb |l6| have shown that intensive monitoring of a microlensed af- 
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terglow lightcurve can be used to reconstruct the parameters of the fireball and 
its environment. The dependence of the afterglow image on frequency offers a 
fingerprint that can be used to identify a microlcnsing event and distinguish it 
from alternative interpretations. It can also be used to constrain the relativistic 
dynamics of the fireball and the properties of its gaseous environment. At the 
highest frequencies, the divergence of the surface brightness near the edge of the 
afterglow image (r = 1 in Figure ||) depends on the thickness of the emitting 
layer behind the relativistic shock front, which is affected by the length scale 
required for particle acceleration and magnetic field amplification behind the 
shock[5l2ltl. 



Ioka & Nakamura 26 considered the more complicated case where the ex- 
plosion is collimated and centered around the viewing axis. More general orien- 
tations that violate circular symmetry need to be considered in the future. 
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Fig. 6. UBVRIJK photometry of GRB 000301C as a function of time in days from 
the GRB trigger (from Garnavich et al. 2000; Gaudi et al. 2001). Data points have 
been offset by the indicated amount for clarity. The dashed line is the best-fit smooth, 
double power-law lightcurve (with no lensing) , while the solid line is the overall best-fit 
microlensing model, where the SBP has been determined from direct inversion. The 
inset shows the i?-band data only. The dotted line is the best-fit microlensing model 
with theoretically calculated SBP, for k = and v > v c . 
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4.1 GRB 000301C 

Garnavich, Loeb, & Stanek Jlf| have reported the possible detection of a mi- 
crolensing magnification feature in the optical-infrared light curve of GRB 000301C 
(see Figure ^). The achromatic transient feature is well fitted by a microlensing 
event of a 0.5Mq lens separated by an Einstein angle from the source center, and 
resembles the prediction of Loeb & Perna |H| for a ring-like source image with 
a narrow fractional width (~ 10%). Alternative interpretations relate the tran- 
sient achromatic brightening to a higher density clump into which the fireball 
propagates j|, or to a refreshment of the decelerating shock either by a shell 
which catches up with it from behind or by continuous energy injection from 
the source |4^| . However, the microlensing model has a smaller number of free 
parameters. If with better data, a future event will show the generic temporal 
and spectral characteristics of a microlensing event, then these alternative inter- 
pretations will be much less viable. A galaxy 2" from GRB 000301C might be 
the host of the stellar lens, but current data provides only an upper-limit on its 
surface brightness at the GRB position. The existence of an intervening galaxy 
increases the probability for microlensing over that of a random line-of-sight. 

Gaudi, Granot, & Loeb [|l7| have shown that direct inversion of the observed 
light curve for GRB 000301C yields a surface brightness profile (SBP) of the 
afterglow image which is strongly limb-brightened, as expected theoretically (see 
Figure Q). 

Obviously, realistic lens systems could be more complicated due to the exter- 
nal shear of the host galaxy or a binary companion. Mao & Loeb |28| calculated 
the magnification light curves in these cases, and found that binary lenses may 
produce multiple peaks of magnification. They also demonstrated that all after- 
glows are likely to show variability at the level of a few percent about a year 
following the explosion, due to stars which are separated by tens of Einstein 
angles from the line-of-sight. 

What is the probability for microlensing? If the lenses are not strongly clus- 
tered so that their cross-sections overlap on the sky, then the probability for 
having an intervening lens star at a projected angular separation 9 from a source 
at a redshift z ~ 2 is ~ 0.3^(9 /9 E ) 2 H!I1||, 

where J?* is the cosmological 
density parameter of stars. The value of .!?* is bounded between the density of 
the luminous stars in galaxies and the total baryonic density as inferred from 
Big Bang nucleosynthesis, 7 x 10 -3 < fl± < 5 x 1CT 2 |]|. Hence, all GRB af- 
terglows should show evidence for events with 9 < 309 e, for which microlensing 
provides a small perturbation to the light curve (This crude estimate ig- 
nores the need to subtract those stars which are embedded in the dense central 
regions of galaxies, where macrolensing dominates and the microlensing optical 
depth is of order unity.) However, only one out of roughly a hundred afterglows 
is expected to be strongly microlensed with an impact parameter smaller than 
the Einstein angle. Indeed, Koopmans & Wambsganss [^7) have found that the 
'a posteriori' probability for the observed microlensing event of GRB 00030 1C 
along a random line-of-sight is between 0.7-2.7% if 20-100% of the dark matter 
is in compact objects. 
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Fig. 7. Fitting GRB 000301C with different SBPs as a function of normalized radius 
(taken from Gaudi et al. 2001). The points are the SBPs determined from direct in- 
version, with la errors (solid) and 3cr errors (dotted). The curves are theoretically 
calculated SBPs for various frequency regimes (see Figure ||). (a) Uniform external 
medium, k = 0. (b) Stellar wind environment, k — 2. The number of degrees of free- 
dom is 92 for the direct inversion points and 89 for the curves. 



Microlensing events are rare but precious. Detailed monitoring of a few strong 
microlensing events among the hundreds of afterglows detected per year by the 
forthcoming Swift satellite, could be used to constrain the environment and the 
dynamics of relativistic GRB fireballs, as well as their magnetic structure and 
particle acceleration process. 
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